Impact of turbulence in long range quantum and classical communications 
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The study of the free-space distribution of quantum correlations is necessary for any future application of 
quantum as classical communication aiming to connect two remote locations. Here we study the propagation of 
a coherent laser beam over 143 Km (between Tenerife and La Palma Islands of the Canary archipelagos). By 
attenuating the beam we also studied the propagation at the single photon level. We investigated the statistic of 
arrival of the incoming photons and the scintillation of the beam. From the analysis of the data, we propose the 
exploitation of turbulence to improve the SNR of the signal. 
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Introduction - The study of the free-space propagation of 
quantum correlations is necessary for any future application 
of quantum communication aiming to connect two remote lo- 
cations. The problem related to the free-space propagation 
is represented by the atmospherical turbulence, that acts as a 
temporal and spatial variation of the air refraction index. A 
turbulent channel acts an increment of the losses on the trans- 
mitted photons due to beam-wandering of the beam-centroid 
or scintillation, increasing the role of the noise fl]-[5]. The un- 
derstanding of the propagation effects induced by turbulence 
at the receiver as well as the temporal statistics of the incom- 
ing photons is crucial to assess the quality of the communi- 
cation and eventually the feasibility of the free-space ground- 
ground and space-ground links j6|-[8). 

In this work we study the propagation of a free space opti- 
cal link (143Km) between Tenerife and La Palma Islands of 
the Canary archipelagos ||9 WTI . The transmitter is located at 
La Palma, on the roof of the Jacobus Kapteyn Telescope (JKT) 
and the receiver is the Optical Ground Station (OGS) at Tener- 
ife. The campaign was performed during the nights between 
17 and 25 September 201 1. 

Optical Setup - The optical setup of the transmitter is 
shown in figure [T] It consists of suitably designed telescope 
whose key component is a singlet aspheric lens of 23 cm di- 
ameter and 220 cm focal lenght at 810nm. The lens diameter 
was chosen to be significantly greater of the estimated Fried 
parameter r$ lfl2l in order to obtain a spot at OGS compared 
to the telescope primary mirror and consequently a greater 
power transfer between the two sites. Our light source is an 
infrared diode at 808nm coupled into single mode fiber with 
output power of about 6mW and suitable attenuators. In or- 
der to facilitate the raw pointing a mechanical XY stage has 
been added (we define the Z direction as the optical axis of the 
system). This stage moves all the 2.5m long telescope in the 
XY plane. All the structure is assembled by three aluminum 
flanges, one hods the lens, one the focal plane and the other is 
attached to the XY back stage. The lens is fixed to a articu- 
lated mount to prevent bending of the structure. 

The IR source has been aligned by means of a dichroic mir- 




FIG. 1: Schematic of the optical setup. 



ror (DM) reflecting the 808nm radiation. The platform carry- 
ing the focusing lens, the collimating lens and the fiber port for 
the IR can be moved by a micro metric XYZ system based on 
stepped motors. In this way the beam can be slightly steared 
by moving the focal spot at the singlet focus. The instanta- 
neous deviation from the initial pointing direction is acquired 
by using a 532nm beacon beam sent from the receiver by us- 
ing a small portable low power laser module directly pointed 
towards La Palma without any optics. The beacon laser is ac- 
quired with a CMOS camera placed on the movable platform 
after the DM trasmission. The centroid of the beacon spot 
on the camera determines the correction on the outcoming IR 
laser by means of an error signal with respect to the reference 
position. The position of the spot at the camera are sampled 
about once a second and averaged for a number of frames and 
this data feeds a control software that calculates the move- 
ment for the fine XY stage in order to compensate slow drift 
in beam direction, which was performed in 1 sec timescale. 

We collected data at the OGS in Tenerife in order to mea- 
sure the received power, the scintillation and analyze the tem- 
poral structure of the signal. We placed in the OGS Coude 
focus a photodiode, a power meter and, when the beam is 
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suitably attenuated with neutral filter, we also collected data 
with a single photon (Excelitas SPCM-AQRH model) detector 
(SPAD). The following data were recorded during the nights 
between 21th and the 24th September 2011. We obtained an 
average attenuation of about 30dB for many times during the 
best run with peaks of 27dB averaged over 2 minutes. The at- 
tenuation is calculated from the fiber and not from the singlet 
lens: the attenuation of the telescope is thus included in the 
measured attenuation. 

Link Analisys - Let's first describe the single photon detec- 
tion acquisition. We performed several measurements by set- 
ting the counting interval T to 0.1ms, 1ms and 10ms. Due 
to turbulence effects, the mean photon number q in a counting 
interval at the receiver should follow a lognormal probability 
distribution [131: 
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where (q) is the average, a 1 = ln(l + SI) and SI = 
is the scintillation index. If the counting interval T is large 
compared with the coherence time of the source and T is short 
compared with the turbulence timescale, the probability of de- 
tecting n photon in each interval follow the Mandel distribu- 
tion: 
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Note that the mean number of detected photon is (n) — 
Tim n P n ~ (l) ■ ^ e report the analysis of the temporal distri- 
bution of an acquisition with 1ms counting interval in figure 
[2jtop). It is possible to observe that, when the average num- 
ber of detected photons, (n), is large (typically larger than 50) 
and the scintillation index bigger than 1, the lognormal and 
Mandel distribution are quite similar. Given the experimental 
scintillation as 2.23±0.01 and the mean value of detected pho- 
ton as 234, we show the counting occurrences together with 
the corresponding lognormal distribution in figure [^bottom). 
For comparison, we also insert the corresponding Mandel dis- 
tribution with the (q) = 234.1 ± 0.1 and a = 349.2 ± 0.2 
parameter obtained from the raw data. We evaluated the simi- 
larity between the experimental data and the lognormal or the 
Mandel curve, defined as S = v^jj where n, and o, are 
respectively the theoretical and experimental occurrence. The 
similarity of the lognormal curve with the data is 0.9959 while 
the Mandel curve has a similarity of 0.9967 showing a clear 
evidence of the statistic transformation. The green curve rep- 
resents the corresponding Poisson distribution with the same 
observed mean value, to compare what it have been obtained if 
the statistic of arrival photon would be purely poissoninan. In 
table|I]we report the data obtained for several different SPAD 
acquisitions. 

As said, we also measure the intensity of received light with 
a fast photodiode by using an intense laser source. In figure 
[3] we plot the temporal distribution of the photodiode volt- 
age of a data set covering 20 s. The intensities are recorded 



Mean= 234.1 ±0.1, std= 349.2±0.2, Scint lndex= 2.23±0.01 
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FIG. 2: SPAD temporal distribution of count occurences and corre- 
sponding lognormal and Mandel curves. We can compare the data 
with the corresponding Poissonian distribution with the same mean 
value (234.1 ± 0.1) that would be obtained without turbolence. 



Mean= 10.6±0.1mV, std= 11.5±0.1mV, Scint lndex= 1.19t0.01 
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FIG. 3: Photodiode temporal distribution intensity occurences and 
corresponding lognormal curve. 



with 50kHz frequency. Also in this case the intensity oc- 
curences follows a lognormal distribution ([TJ as shown from 
the lognormal curve with similarity of 0.9896. In this case 
the scintillation index evaluated from the experimental data is 

SI = fj^r = 1.19 ±0.01. 

Improving the SNR - For both quantum and classical com- 
munication, it is of paramount importance achieving an high 
signal to noise ratio (SNR). If a qubit state l^) encoded in the 
photon polarization must be sent between two remote loca- 
tion, it is possible to determine the effect of (white) noise on 
the polarization fidelitv lfT4l . Let's measure the SNR in dB, 
namely SNR= 10 log 10 ^r-, where N n the average amount of 
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FIG. 4: SPAD power spectrum and cumulative power spectrum. 
Frequency bound (red line): the frequencies below 51.1724 Hz con- 
tribute to the 95% of the scintillation index. 



noise (coming from dark detections or background radiation) 
and N s is the average counting signal. It is easy to show that 
the fidelity depends on the SNR as: 
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In order to improve the SNR for the transmission of sin- 
gle photons in a long distance free-space link as the present 
one which use a lm optical receiver, out of our findings we 
envisage the exploitation of the following procedure. With a 
given frequency (slower than the single photon transmission 
rate) the free-space channel is probed by means of a classical 
signal that gives the information of the instantaneous trans- 
mission of the channel. Only if the transmission is above a 
given threshold the single photon signal is acquired. It is cru- 
cial for the protocol to be efficient to correctly identify the 
"probing" frequency and the threshold to be used. This tech- 
nique can be also used in the classical case, for instance in the 
on-off keying. 

We report in figure|4]the frequency spectrum and the cumu- 
lative power spectrum of the data plotted in figure [2] The nor- 
malized plot of the power spectrum is obtained by normalizing 
the intensities by the average I' — 1/(1). The power spectrum 
is related to the scintillation index as follow. We write the set 
of (normalized) acquisitions as II with k = 0, • • • , N— 1 and 
N = 20s/20/is = 10 6 the number of intensity acquisitions 
over the 20 seconds. The Fourier components are given by 
/„ = Ylk l'k ujnk w i tn w = e^^r . By Parseval's theorem it 
is easy to show that SI — ~Y^1=\ | -^n | 2 , namely it is the 
cumulative power spectrum without the zero frequency (Iq) 
component. We can notice that the frequencies contributing 
to the scintillation (up to 95%) are within (almost) 50Hz. For 
frequencies above around 500Hz, the spectrum becomes flat, 
indicating that at this frequency the random noise is domi- 
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FIG. 5: Duration (in ms) of events with over-threshold counting: in 
the different plots we considered a threshold of 1 , 2, 4 and 6 dB above 
the average. 



nant. The typical fluctuations of the transmission channel due 
to turbulence are within 100Hz (see table The frequency 
analysis of the temporal scintillation indicates that the prob- 
ing frequency doesn't need to be higher than lKHz. 

In order to obtain a further evidence, we analysed the fea- 
tures of the counts above a given threshold of the signal re- 
ported in fig. [2] By considering a threshold of 1, 2, 4 and 6 dB 
above the average we considered the duration (in ms) of events 
with over-threshold counting. The results are shown in figure 
[5] The probability of obtaining an event above a given thresh- 
old q can be predicted from the lognormal distribution! 15 1 
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where erf(a;) is the Gaussian error function erf(x) = 
J* e~ t dt. Acquiring the single photon channel only if 
the "probed" transmission is above a given threshold implies 
an increase of the average photon counts in each time-slot. 
It is possible to show that, by considering only the events in 
which the transmission satisfy T > Tq, the new mean value 

{n)thres IS 



(^7 thres 



{II 



1 - erf [ 



In - 



> 1 



(5) 



1 - erf [- 



/2cr 2 



4 




Threshold (dB) 

FIG. 6: SNR and the percentage of the overall counts that will be 
detected in function of the threshold selection. 



Clearly, this threshold selection increase the SNR, but at the 
same time decrease the overall counts in given time. In figure 
|6]we show the increase (in dB) of the SNR and the percentage 
of the overall counts that will be detected. In cases of strong 
turbulence and high noise, this technique could help the qubit 
transmission by "exploiting turbulence", namely considering 
only the particular moments in which the turbulence increase 
the channel transmission. 

In conclusion, the statistic of arrival of single photon over 
free-space 143 Km optical link has been analyzed demonstrat- 
ing the transformation from Poissonian to Lognormal distribu- 
tion thus expanding this investigation for more than an order 
of magnitude in length with respect to previous results lfl3l . 
The evidence of consecutive subintervals of low losses allows 
to envisage the exploitation of turbulence as an SNR improve- 
ment technique. 
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1.5485 
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132.10 
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Low 


37.86 


48.42 


200 
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1.6355 


35 


Low 


20.83 


24.18 


200 
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1.3479 


35 
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2.862 


3.795 


65 
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1.7578 


384 


Low 


5.267 


7.519 


65 


1 


2.0383 
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TABLE I: Data obtained for different single photon acquisition com- 
pared to the background (first line). For each acquisition we report 
the total duration of the acquisition (Time), the temporal windows 
defining the counting interval (Window), the mean number of counts 
in the counting interval (Mean) and its the standard deviation (Std). 
We also report the scintillation index (SI) and the frequency bound 
such that all the frequencies below the bound contribute up to 95% of 
the scintillation (Bound). With High (Low) we indicate acquisition 
with high (low) mean photon number detected during Is. We notice 
that for the last two data sets the bound is higher due to the low signal 
compared to the background (having flat frequency spectrum). 
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